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Abstract 

We discuss the properties of fluctuations of the electric charge in the vicinity of 
the chiral crossover transition within effective chiral models at finite temperature 
and vanishing net baryon density. The calculation includes non-perturbative 
dynamics implemented within the functional rcnormalization group approach. 
We study the temperature dependence of the electric charge susceptibilities in 
the linear sigma model and explore the role of quantum statistics. Within the 
Polyakov loop extended quark-meson model, we study the influence of the cou- 
pling of quarks to mesons and to an effective gluon field on charge fluctuations. 
We find a clear signal for the chiral crossover transition in the fluctuations of 
the electric charge. Accordingly, we stress the role of higher order cumulants 
as probes of criticality related to the restoration of chiral symmetry and decon- 
fincment. 

Keywords: QCD phase diagram; Heavy ion collisions; Chiral phase transition; 
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1. Introduction 

Thermodynamics of strongly interacting matter has been explored numeri- 
cally within Lattice Quantum Chromodynamics (LQCD) [H, 0, E| H| as well as 
experimentally in heavy- ion collisions |5| . The fundamental goal of these stud- 
ies is to reveal the structure of the QCD phase diagram and to investigate the 
critical properties of QCD. 

LQCD results show that at finite temperatures QCD exhibits restoration 
of chiral symmetry and deconfinement in a crossover transition. The LQCD 
equation of state indicates a clear separation between the confined hadronic 
phase and the deconfined quark-gluon plasma. The temperature separation of 
the phases of QCD is very transparent when one considers fluctuations of con- 
served charges @, 0, 0, 0, 0, [Io[ • At vanishing chemical potential the variance 
of both the net baryon number and the electric charge exhibits a rapid change 
in a narrow temperature interval. This behavior is attributed to the deconfine- 
ment of quarks [Til IT2I [l3| . On the other hand, the fourth order cumulants of 
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these charges show peak-like structures, whereas the sixth order cumulants are 
negative in a narrow temperature interval [(| . 

These properties of the fluctuations observed in LQCD at finite temperature 
and for small masses of the up- and down-quarks can be attributed to the 
critical dynamics of the explicitly broken chiral symmetry, which is characterized 
by the 0(4) universality class of QCD As noted in [13, El, this could have 

interesting phcnomcnological implications: it has been suggested, that the sixth 
and eight order cumulants of the net baryon number and electric charge may 
allow an experimental verification of the chiral crossover transition in heavy-ion 



collisions 15 



The critical behavior of strongly interacting matter related with chiral dy- 
namics should be common to all models which exhibit the underlying chiral 
symmetry of QCD and the same patterns of spontaneous chiral symmetry break- 
ing at finite temperature. Thus, such effective models have been used to study 
the thermodynamics near the chiral phase transition and to explore observablcs 
which are sensitive to the critical behavior [l6, 17, 18, 19|, 2(J 21, 22, 23 1. In this 
context, models which include the coupling of quarks to an effective gluon field, 
such as the Polyakov loop extended Nambu-Jona-Lasinio (PNJL) [18j, [2l| and 
the quark-meson (PQM) [24[ models, are very useful. Both models reproduce 
the essential features of the QCD thermodynamics already in the mean-field 
approximation [l|| [13, EH, Il9l.l2l|. However, to correctly account for the crit- 
ical dynamics, it is necessary to go beyond the mean-field approximation and 
include non-perturbative fluctuations. This can be achieved e.g . by e mployin g 



methods based on the functional renormalization group (FRG) [25l. l26t 1271 l28j. 
The FRG has been applied to study the fluctuations of the net baryon number 
near the chiral phase transition at finite and at vanishing baryon chemical po- 
tential [l4, 29|, 11, 13, 3(J. However, so far such studies are not available for the 
electric charge fluctuations. 

In this work, we study fluctuations of the electric charge near the chiral phase 
transition at finite temperature and at vanishing chemical potential within the 
FRG approach to the linear sigma and PQM models. We explore the tem- 
perature dependence of these fluctuations and the role of quantum statistics. 
Furthermore, the influence of quark- and gluon-dynamics on the properties of 
the electric charge fluctuations near the chiral crossover transition, including the 
effect of mesonic fluctuations, is discussed. We find a strong sensitivity of higher 
order cumulants of the electric charge to chiral symmetry breaking and decon- 
finemenet. We also show, that these cumulants exhibit a characteristic 0(4) 
structure, similar to that found in the fluctuations of the net baryon charge. In 
particular, the cumulants turn negative in the vicinity of the chiral transition 
also at vanishing net charge. 

In the subsequent section we introduce the non-perturbative thermodynamic 
potential in the FRG approach to the PQM model formulated at finite electric- 
charge chemical potential. In section 3 we present our results on fluctuations of 
the electric charge, including higher cumulants, in the linear sigma and PQM 
models. Finally, in section 4 we present our conclusions. 
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2. The thermodynamic potential in the Polyakov-quark-meson model 



The quark-meson model is an effective realization of the low-energy sec- 
tor of QCD, which incorporates chiral symmetry and exhibits a global SU(N C ) 
color symmetry. Hence, it docs not describe quark confinement. Neverthe- 
less, by introducing a coupling of the quarks to a uniform temporal color gauge 
field, represented by the Polyakov loop, one obtains an equation of state with 
properties that closely resemble those of QCD matter, including a change of 
effective degrees of freedom with increasin g te mperature, thus mimicking the 
confinement-deconfinement transition 17, ljl 21, 24|. This is sometimes re- 
ferred to a statistical confinement. 

A formulation of the thermodynamics of the PQM model, which remains 
valid near the chiral phase transition, requires the use of non-perturbative meth- 
ods. In an analysis of electric charge fluctuations, it is of particular importance 
to account for mesonic fluctuations, including specifically interacting pions as 
carriers of electric charge. We employ a method based on the functional renor- 
malization group (FRG) to compute the thermodynamic potential in the PQM 
model. This method involves an infrared regularization of the fluctuations at 
a slidin g m omentum scale k, resulting in a scale-dependent effective action 
[H, I2II |27|. 28 1 . We treat the Polyakov loop as a background field, which is 



introduced self-consistently on the mean-field level while the quark and meson 
fields, fluctuations are accounted for by solving the FRG flow equations. 

We follow the procedure used in Ref. [2^] in the formulation of the flow equa- 
tion for the scale-dependent grand canonical potential density, Q k = TTk/V, for 
the quark and meson subsystems at finite temperature and for a non-vanishing 
electric charge chemical potential. The thermodynamic potential is obtained by 
solving the flow equation 
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Here n B (E;T) is the bosonic distribution function 
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n B {E-T) = 
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= e^fXQ is the charge pion chemical potential and = 1 the charge of a n + . 
The pion and sigma energies are given by 



E n = \lk 2 + n' k , E a = Jk 2 + n' k + 2 P n" 
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where the primes denote derivatives with respect to p = (a 2 + t? 2 )/2 of Q = 
Q + ca. The fermion distribution functions N(£,£*;T, fj, f ) and N(£,£*;T, fj, f ), 

N(£,£ ;T,fj, f ) = - - UemEq _ M) - u * e p( Eq - M ) + e 3/3(B,- M/ ) ' ( 2 ) 

N(£,£*;T,v f ) = N(£*,£;T,-» f ), (3) 

are modified because of the coupling to the gluon field. Finally, the quark energy 

reads 

E q = ^WTWp (4) 
and the quark chemical potentials arc defined by 

1 1 

Mm = g MB + e nMQi Md = g Ms + e d MQ (5) 

with e u = 2/3 and e<2 = —1/3. 

The flow equation |T]) is solved numerically with the ultraviolet cutoff A = 1.2 
GeV using the polynomial method described in Ref. [29[. In this scheme, the 
stationarity condition 



da 



da 



(6) 



is implemented in the flow equation. The initial conditions for the flow are 
chosen to reproduce the following vacuum properties: the physical pion mass 
m n = 138 MeV, the pion decay constant / ff = 93 MeV, the sigma mass m a = 
600 MeV, and the constituent quark mass m q = 300 MeV at the scale k — > 0. 
The symmetry breaking term, c = m^/^, is treated as an external field which 
does not flow. The flow of the Yukawa coupling g is neglected because it is not 
expected to be significant for the present studies (see e.g. Refs. [3ll|). 

By solving Eq. {l} one obtains the thermodynamic potential rtk->o(£, £*;T, /x) 
as a function of the Polyakov loop variables £ and I* . The full thermodynamic 
potential f2(£, I*; T, /j,) in the PQM model, including quark, meson and gluon 
degrees of freedom, is obtained by adding the effective gluon potential U(£,£*), 

Q{£, t;T, n) = n k ^ (£, t;T, /x) + U{£, £*). (7) 

At a given temperature and chemical potential, the Polyakov loop variables, I 
and £*, are then determined by the stationarity conditions: 

J^,r ; T,/i) = o, A^,r ; i» = o. (8) 

The thermodynamic potential (J7J does not contain contributions of thermal 
modes with momenta larger than the cutoff A. In order to obtain the correct 
high-temperature behavior of the thermodynamic functions, we supplement the 
FRG potential with the contribution of the high-momentum states. A procedure 
for implementing this was proposed in Ref. [321 ] for the QM model and extended 
to the PQM model in Ref. [29(, by including the flow of quarks interacting with 
the Polyakov loop for momenta k > A. 
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3. Electric charge density fluctuations 



The thermodynamic potential obtained by solving the flow equation ((T|) cor- 
rectly reproduces the critical 0(4) scaling near the chiral phase transition. Thus, 
using this potential one can compute the temperature dependence of the electric 
charge fluctuations and their higher cumulants. 

The fluctuations of the electric charge are characterized by the generalized 
susceptibilities, 

Q(T) = d"b(J> Q )/r*] 

The first cumulant, x? = fig/T 3 , is given by the electric charge density ng = 
Nq/V, while the second cumulant 

X% = ^((SN Q f) (10) 

with 5Nq = Nq — (Nq), is proportional to the variance of the electric charge. 
Furthermore, the fourth and sixth order cumulants can be expressed through 
moments of SNq 

X? = T^[((SN Q ) 4 )-3{(SN Q fr], (11) 



X? = TFF3 [((5N Q ) e )-15((5N Q r)((SN Q f)~10((SN Q ) 3 f (12) 
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1 
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Based on LQCD results [33|, one expects that in the region of the crossover 
transition, the fluctuations of the electric charge should reflect the critical scal- 
ing of the underlying 0(4) symmetry. For the 0(4) universality class, the sin- 
gular part of thermodynamic pressure, at small values of the baryon chemical 
potential, scales as 

pcx(T/T c -l + ^!) 2 - Q , (13) 

where T c is the temperature of the second-order phase transition at fi = in 
the chiral limit. The exponent a « —0.21 is the critical exponent of the specific 
heat in the 0(4) universality class in three-dimensions. Therefore, at vanishing 
chemical potential, the cumulants Xn 01 the net baryon number diverge at the 
transition temperature for all even n > 6, 

X f„cx(T/r c -l) 2 -"- Q . (14) 

Fluctuations of the electric charge are related to the net baryon and isospin 
fluctuations by 



y q = — 

An 2?l 



Xn Xn ~t~ ^ ^ 



9 n (p/3 4 ) 



(15) 
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Figure 1: (Color online) The second \2 anc ' * nc fourth \a or der cumulant of the net electric 
charge fluctuations in the linear sigma model obtained within the functional renormalisation 
group approach (solid lines). The dotted line is the contribution of an ideal pion gas in the 
Boltzmann approximation. The dash-dotted (l ma x = 1) and dashed (£ ma;c = 2) lines show the 
first and the second terms in the series 1171 . respectively. The temperature axis is normalized 
by T pc = 180.5 MeV which is the pseudo-critical temperature in the model within the FRG 
approach. 



The isovector fluctuations are regular at the chiral phase transition in an isospin 
symmetric system. Thus, the xli remain finite at the phase transition. The last 
term, which correlates isospin with baryon number, can diverge for sufficiently 
large n. However, these terms give subleading contributions to criticality of \n- 
Thus, the leading singular part of the electric-charge cumulants Xn ^ solely 
determined by Xn- Consequently, at small quark masses, the contributions 
from the singular part of the thermodynamic pressure to the net electric charge 
fluctuations should be similar to that found for the net baryon fluctuations. 

The expected residual 0(4) scaling in Xn & t finite quark masses is indeed 
observed in the PQM model in the FRG approach. However, since the regular 
contributions to the electric charge and to the baryon number fluctuations are 
different, it is not a priori clear that this will be the case also for cumulants Xn 
of the electric charge. 

In the following, we compute the temperature dependence of Xn m ^ ne PQM 
model with the thermodynamic potential obtained from Eq. ([7]) . We also discuss 
the contributions of individual modes to the properties of Xn near the chiral 
transition. 

3.1. Fluctuations in an interacting pion gas 

We first consider the contribution of mesonic degrees of freedom to the fluc- 
tuations of the net electric charge. Clearly, in such studies pionic interactions 
should be included. A suitable model for exploring the meson contribution is 
the mesonic linear sigma model. The thermodynamic potential for this model 
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Figure 2: (Color online) The cumulants of the electric charge fluctuations in the Polyakov loop 
extended quark-meson model calculated within the functional rcnormalisation group approach 
(solid lines). The results are also shown for only charged quarks = and only charged 
pions ef = (see text). The hadron resonance gas contribution is shown as dotted lines. 



is obtained directly from the flow equation ([7]) in the limit of N c —> 0. The 
fluctuations of the electric charge are then due only to the charged pion degrees 
of freedom. 

In the mean-field approximation to the linear sigma model, fluctuations of 



the electric charge are absent. They can be included perturbatively [34|, |35 



36j or non-perturbatively within the FRG approach, where fluctuations due 



to interacting pions are included [11|, [12|, |29|, [3Jj . Here we employ the FRG 



approach to account for meson fluctuations. 

The second and the fourth order cumulants obtained in the FRG approach 
are shown in Fig. [1] Although the chiral symmetry is explicitly broken by the 
vacuum pion mass, a clear signature of the chiral crossover transition is seen 
in the susceptibilities. In particular, we note the non-monotonic temperature 
dependence of the fourth order cumulant, resulting in a peaked structure near 
the chiral transition. This behavior is due mostly to the temperature depen- 
dence of the pion mass, which up to the chiral transition remains close to its 
vacuum value, while above the transition grows strongly with temperature. 
Consequently, above the pseudo-critical temperature (T pc ), all cumulants of the 
net-charge fluctuations decrease with temperature. Note that this effect is not 
a consequence of the cut off, which suppresses high momentum modes in the 
model. Due to the large value of A = 1.2 GeV, this suppression sets in only at 
much higher temperatures. 

It is interesting to asses at which temperature pionic interactions start influ- 
encing the charge fluctuations. This can be estimated by comparing the FRG 
results for the X2 an d Xa w ^h those obtained for an ideal gas of charged pions. 
The pressure of the latter is given by 

= 2^7* E T* K2{Wm) ^ + ^ ' (16) 
1=1 

where Ki (x) is a modified Besscl function and [i q is the electric chemical poten- 
tial. In an ideal pion gas, the odd-order cumulants of the electric charge vanish 
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at \x q = 0, while even cumulants are given by 

2 00 

^ = iliE r ^( ; H. (17) 

2 = 1 

This can be considered as an approximation to the regular part of the fluctua- 
tions of the electric charge. 

In Fig. [1] the contributions of the regular part to the second and fourth 
cumulants are shown. At low temperatures the FRG and ideal pion gas results 
agree. This indicates that the FRG method correctly accounts for the pionic 
contribution to the thermodynamics. 

The effect of quantum statistics can be studied by truncating the series in 
Eq. (fT7]) at a given order I = l max . This is illustrated in Fig. Q] where we 
confront the l max — 1 and 2 truncations for the second and the fourth order cu- 
mulants of the free pion gas with the full (l ma x = 00) result. Clearly, the leading 
contribution, which corresponds to Boltzmann statistics, deviates considerably 
from the full quantum statistics results, except at very low temperatures. The 
higher-order terms in the series contribute effectively as multiply charged par- 
ticles. Hence, their contribution increases with the order of the cumulant. This 
expectation is confirmed by Fig. [TJ which shows that the deviation is stronger 
for Xa than for X2 ■ Thus, the effect of quantum statistics cannot be neglected 
in the calculation of fluctuations of the electric charge. This is a reflection of 
the fact that multi-charged particles can yield a dominant contribution to the 
fluctuations of conserved charges [f|. In QCD this is the case for fluctuations of 
the electric charge and strangeness, where double- and triple-charged baryons 
play an important role. 

3.2. Electric charge cumulants in the PQM model 

The dynamics of the PQM model and its particle content is more relevant for 
QCD than the mesonic linear sigma model described above. The PQM model 
contains chargcd-quarks as dynamical degrees of freedom and exhibits statistical 
confinement through the coupling of the quarks to the effective gluon field. 

In Fig. [5] we show the first three non-trivial cumulants X2 'j Xa an d Xe °f 
the electric charge for vanishing baryon- and electric charge chemical potential 
obtained from Eq. ([7]). We also show separately the contributions of quarks or 
pions. This separation is made by setting the electric charge of pions e T and 
then that of quarks ef to zero, respectively. 

The second-order cumulant X2 snows a rapid increase near the pseudo- 
critical temperature, interpolating between the pion and quark contributions. 
Since X2 1S n °t influenced by the singular part of the thermodynamic pressure 
(fP3|) . the behavior seen in Fig. [2] is due to the "statistical confinement" of the 
PQM model. The rapid increase with temperature is a consequence of the rapid 
unleashing of single- and double-quarks states at the crossover temperature. 
The pion contribution, on the other hand, is suppressed above T pc , owing to the 
increasing thermal mass. 
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Figure 3: (Color online) The kurtosis (X4/X2) ant ^ Omi /X? ) ratio calculated in the 
Polyakov loop extended quark-meson model. The arrow shows the corresponding Stefan- 
Boltzmann limits. The hadron resonance gas model results are indicated as dottcd-lincs. 



The fourth- and sixth-order cumulants exhibit a peak, which increases in 
strength with the order of the cumulant. This can be understood in terms 
of the critical dynamics of the chiral transition; with increasing order of the 
cumulant, the contribution of the singular part of the pressure (|13|) becomes 
more and more dominant. The negative structure of the Xe near T pc is similar 
to that observed for the sixth order cumulant of net baryon number fluctuations 
and is due to the particular form of the 0(4) scaling function 14|. The peak 
in X4 appears as a sum of quark and pion contributions, while the peak in x& 
is dominantly due to pions. The negative region of Xe near T pc is due to the 
quark contribution. 

In Fig. [2]we also compare the PQM results with that of the hadron resonance 
gas (HRG) model, which includes the contributions of all charged hadrons and 
resonances. The HRG model reproduces the thermodynamics of LQCD in the 
hadronic phase up to T ~ 0.9T pc . 

There is good agreement between the PQM and HRG model results at low 
temperatures, where in both cases pions are the dominant degrees of freedom. 
For higher temperatures the HRG overshoots the PQM model results. Thus, 
the PQM model does not yield a quantitative description of the LQCD results 
on thermodynamics. Nevertheless, this model can provide useful insights into 
the critical dynamics of universal quantities near the chiral phase transition. 

In order to reduce the contribution of the non-singular part to the fluctu- 
ations and to focus on the critical behavior, ratios of Xn t° the second-order 
cumulant X2 j which is not influenced by the critical chiral dynamics at /i = 
[H, have been studied. In Fig. [3] we show the temperature dependence of the 
kurtosis k = X4 /X2 an d of the ratio Xe IX2 near the chiral crossover transi- 
tion. At low temperatures, k — > 1, as expected for a non-interacting Boltzmann 
gas. For higher temperatures, the kurtosis increases owing to the contribution 



of double-charged baryons and to higher order quantum corrections to the pion 
distribution. Near T pc and above the quarks dominate, because of the increas- 
ing pion mass, resulting in k decreasing towards the Stefan-Boltzmann limit of 
free quark gas. The ratio xi Ixi exhibits a much stronger peak, which diverges 
at the critical temperature in the chiral limit. In analogy with the net baryon 
fluctuations [14j, [l5|, this ratio also develops a negative region near T pc , owing 
to the chiral critical dynamics encoded in the 0(4) scaling function. 

As shown in Fig. [3J in the PQM model both ratios X4/X2 an d Xs '/X2 
overshoot the HRG model predictions. This indicates that in the HRG model, 
the contribution of multi-charged states to fluctuations is reduced relative to 
that of singly-charged particles more than in the PQM model. 



4. Summary and conclusions 

Wc have computed the electric charge fluctuations in the vicinity of the chiral 
crossover transition within effective chiral models. Our calculations include 
non-perturbative dynamics implemented within the functional renormalization 
group approach at finite temperature and at vanishing baryon density. We 
have performed our studies in the linear sigma model and in the Polyakov-loop 
extended quark-meson model. 

We find a strong sensitivity of the fluctuations of the electric charge to the 
chiral transition. We have also explored the influence of multi-charged states 
and quantum statistics on the temperature dependence of cumulants of the 
electric charge. 

Our results show that at physical pion and quark masses, the higher order 
cumulants of the electric charge fluctuations exhibit properties which can be 
linked to the chiral critical dynamics expected in the 0(4) universality class. 
The qualitative structure of fluctuations in the PQM model is similar to that 
obtained in lattice QCD. This confirms the phenomenological importance of 
electric charge fluctuations, in particular the sixth- and higher-order cumulants, 
as signatures for the chiral crossover transition in heavy-ion collisions. 
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